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. . . We discuss the new data for the production of the tp' meson in pA collisions at 450 GeV at CERN- 

I SPS (of the NA50-coUaboration) [1]. We extract from the CERN data a{ip'N) ^ 8 mh under the 

, assumption that the i/;' is produced as a result of the space-time evolution of a point-like cc pair 

which expands with time to the full size of the charmonium state. In the analysis we assume the 
existence of a relationship between the distribution of color in a hadron and the cross section of its 
interaction with a nucleon. However, our result is rather sensitive to the pattern of the expansion 
of the wave packet and significantly larger values of a{tp'N) are not ruled out by the data. We show 
that recent CERN data confirm the suggestion of ref. [2] that color fluctuations of the strengths in 
charmonium-nucleon interaction are the major source of suppression of the J/tp yield as observed 
at CERN in both pA and AA collisions. 
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I. INTRODUCTION 

\o ' 

, The aim of the present paper is to show that the theoretical description of the production and annihilation yield of 
' particles with hidden charm, which accounts for color fluctuations in charmonium-nucleon interactions, agrees with 
[ new data of the NA50-collaboration [1] who observed quite different J/tp and ■(/''-nucleon cross sections. Analogous 
■ results have been found at negative Feynman x (xp) in proton- nucleus collisions at 800 GeV at Fermilab by the 
' E866 collaboration [3]: Both experiments found that the charmonium-nucleon cross section is smaller in the target 
I fragmentation region than at midrapidity. This is in good agreement with models that assume that the tp' is produced 
due to a space-time evolution of colorless, point-like cc pairs, which expands with time to its full size and that there 
^ , exists a relationship between the distribution of color in a hadron and the cross section of its interaction with a 
nucleon. The agreement of this scenario with the negative x p Fermilab data [3] with was first demonstrated in ref. [2] . 

The major new effect is that QCD predicts different interaction cross sections for the interaction of different 
charmonium states (J/V'j "0' ^-nd x) with usual hadrons due to their different spatial size if the cross section is 
proportional to the distribution of color in the projectile. This relationship is proved in pQCD (perturbative QCD) [4] 
as equivalent of the QCD evolution equations and it is a plausible suggestion in the nonperturbative QCD regime, 
which assumes smooth matching with pQCD predictions. If this is so, the ip'-rmcleon cross section is expected to be 
on the order of the interaction cross section of nucleons with mesons built of light quarks, e.g. the (f> or the K. 

It is a well known fact that spatially larger hadrons have larger interaction cross sections. For example, the 
inelastic J / ip-micleon cross section was found to be around 3.5 mb at SLAG [5] while the inelastic 7r-nucleon cross 
section at these energies is approximately 20 mb. From charmonium models, e.g. in ref. [6,7] it is known that the 
different charmonium states (J/'0, X ^-nd ip') have different spatial sizes. The average relative distances a/ (r^) with 
(r2) = J ^^)^(r) ■ r^(Pr of the QQ-pairs are a/PIj/^ = 0.38 fm, ^/Jj^^ = 0.57 fm and a/P)^, = 0.76 fm for the 
wave functions 5'(r) of ref. [6]. 

The relation between the spatial size of a color dipole and the interaction cross section was found in pQCD [4] 
for a spatially small dipole as another form of the QCD evolution equations. It will be shown in sect. HE that the 
charmonium-nucleon cross section predicted from pQCD is significantly smaller than the measured ones. We concluded 
in ref. [2] that these cross sections are dominated by non-perturbative contributions. A plausible parametrization 
within the constituent quark model for this soft contribution is reviewed in section II D. 

In ref. [2] it was shown that the production of J/ip^s in pA collisions can be understood if one takes into account the 
production and the subsequent decays of higher mass charmonium state resonances (x,'0') info J/ip's. Those higher 
resonances have larger cross sections for the scattering off nucleons due to their larger spatial size. This leads to a 
significant increase of the effective absorption of J/ip's as compared to the propagation of pure J/V>states. In ref. [8] 
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it was shown also that the production of J/tp^s in AA colhsions is additionally suppressed by the final state interaction 
of charmonium states with newly produced hadrons like tt's, p's and so on. These interactions of charmonium states 
with light hadrons are predominantly soft in the SPS-regime, because they take place at lower energies than the 
collisions with the initial nucleons. 



II. MODEL DESCRIPTION 

A. Semiclassical Glauber approximation 

The suppression factor S for the charmonium production in minimum bias pA collisions can be evaluated within 
the semiclassical approximation (cf. [9]) as 

= A^-ai^N^X) =jl^'Bdzp{B,z)- exp <t{XN)p{B, z')Az^ . (1) 

p{B,z) is the local nuclear ground state density (we used the standard parametrization from [10]). a{XN) is the 
interaction cross section of the charmonium state X with a nucleon. We want to draw attention to the fact that this 
cross section changes with time due to the space-time evolution of color fluctuations. Therefore, it is necessary to 
keep a under the integral. In principle one should deal with the expansion effects on the amplitude level. However, 
numerically the difference is small on the scale of the uncertainties of the modeling of the effect. 

We make the standard assumption that the production of cc pairs is a hard process and that the QCD factorization 
theorem is applicable, similar to the Drell-Yan process. In practice this means for Drell Yan pair production that 
5=1, because a{pA = A ■ a{pp l^lr). Deviations from 5 = 1 for the Drell-Yan process can result 

from nuclear effects on the parton distribution only, which are neglected, because in the kinematics of SPS they are a 
small correction. However, in contrast to the dilepton pair produced in the Drell-Yan process, the charmonium states 
have strong final state interactions. They can be split into open charmed hadrons by colliding with other hadrons. 
Therefore, < 5 < 1 for the charmonium states. 

The S = 1 assumption on the level of the charm quark production leads to a restriction on the range of applicability 
of the discussed model as the gluon distribution is modified at small x due to the shadowing and antishadowing effects. 
At the very least the requirement is that for very high energies 

' 2xo ■ xpniN 

where xo ^ 0.03 is in the region where gluon shadowing sets in. a;p is the Bjorken x of the projectile. For Ucm. = 
the restriction corresponds to xq = TOjy^/2y^, because xp = xt with xt is the Bjorken x of the target. In addition, 
as soon as xa of the nuclear gluon becomes small enough the cc pair is produced at a distance ^ ll{2xA'mi^) from 
the interaction point and one has to take into account the evolution of the pair before it reaches the interaction point. 
Hence the picture of the formation of cc states in the proton fragmentation region is qualitatively different from that 
in the nuclear fragmentation region. 

The suppression factor S of J/'^'s produced in the nuclear medium is calculated as: 

5 = 0.6 • (0.92 • S-^l''' + 0.08 • 5*') + 0.4 • . (3) 

S'^ are the respective suppression factors of the different pure charmonium states X in nuclear matter. Eq. (3) 
accounts for the decay of higher resonances after they left the target nucleus into J/ip's. The fractions of J/'if's that 
are produced in the decays of higher resonances in eq. (3) are taken from ref. [11]. However, in ref. [11] it is assumed 
that the different charmonium states interact with nucleons with the same cross section, which is in disagreement 
with the data from the refs. [1,3]. 

In line with the above discussion we want to stress that Eqs. (1) and (3) are applicable at CERN energies for central 
and negative rapidities, but have to be modified if applied already at jjcrn. ^ at RHIC or higher energies, because 
at higher energies charmonium states can be produced outside of the nucleus and the cc pairs propagate through the 
whole nucleus without forming a hadron. 



Data are often presented in the form: 



The relation between S and a is 



CpA = o-pp ■ A" (4) 
S = A"-'^ . (5) 
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B. Color Fluctuations in Charmonium Rescattering 



The first evaluations of atotiJ/ip — N) have been obtained by applying the Vector Dominance Model (VDM) to 
J/ 4' photoproduction data. This leads to dj/^jq ^ 1 mb for Einc ~ 20 GeV. However, the application of the 
VDM leads to a paradox [12] - one obtains crtot{^' — N) 0.7 • atot{J/4' — although, on the other hand, 
r*' w 2rj/^ in charmonium models like in the refs. [7,6]. This clearly indicates that the charmonium states produced 
in photoproduction are in a smaller than average configuration. Therefore, the VDM significantly underestimates 
atotiJ/4' - N) and atoti^' - N) [12]. 

Note that the Generalized VDM, which assumes that the dominant process is photoproduction of spatially small 
cc pairs and accounts for the space time evolution of this pair, predicts significantly larger atoti'^lj' — N) [13]. 

Indeed, the ^-dependence of the J/ip production studied at SLAG at Einc ~ 20 GeV exhibits a significant absorption 
effect [5] corresponding to aabs {J /i' — N) = 3.5 ± 0.8 mb. It was demonstrated in [14] that, in the kinematic region at 
SLAG, the effects due to the space-time evolution of the J/ip are still small for the formation of J ftp's and lead only 
to a small increase of the value of <Jabs{J/fp — N). So, in contrast to the findings at higher energies, at intermediate 
energies this process measures the genuine J/iJj — N interaction cross section at energies of ^ 15-20 GeV [14]. However, 
the dynamical effect of the production of charmonium states in squeezed configurations is still there. We account for 
this eff'ect as due to the propagation of the J/ip and ip' system. Experimental evidence for this expansion was found 
e.g. in the photoproduction of charmonium states. This will be discussed in sect. II C. 

In the semiclassical Glauber approximation, we take into account these color fluctuations in an effective way as 
described in ref. [14]. We assume that a charmonium state X is produced at z as small cc configuration, then it 
evolves - during the formation time tf resp. while it passes the formation length Ic - to its full size. Please note that 
there is up to now no theoretical or experimental proof for the assumption that charmonium states are produced in 
point-like configurations as predicted in PQGD, a way to test this experimentally was suggested recently in ref. [15]. 
Therefore, if the formation length of the charmonium states. If, becomes larger than the average internucleon distance 
{If > rjvAT ~ 1-8 fm), one has to take into account the evolution of the cross sections with the distance from the 
production point z' — z [14]. 

— Z 

a{z, z')x = (7{z) -\ {ax - a{z)) for z' - z <tf 

h 

a{z, z')x — crx otherwise. (6) 

The formation length of the J ftp is given by the energy denominator I f w — -i "^^ -i — , where p is the momentum of 

the J/tp in the rest frame of the target. With p = 30 GeV, the momentum oi a, J/tp produced at midrapidity at SPS 
energies {Eiab = 200 AGeV), this yields Z/ « 3 fm, i.e., a proper formation time of r/ = 0.3 fm. 

As formation time of the ij/ in its rest system we use the radius given by nonrelativistic charmonium models, e.g. 
see the refs. [7,6]. This radius is r = 0.45 fm for the ip'. Please note that this radius differs from the variable r in 
charmonium models, which rather denotes the diameter of positronium like states. Within the formation time the 
cross section increases linearly with the distance from the production point [14,2]. 

We chose a larger value for the cross section of the ijj' N interaction than for the J/ tp because the radius of tp' is 
a factor two larger and its radius is larger than the formation time of the J/ip. A larger value of for the tp' is 
supported also by the extraction of the formation time of the J/tp [16]. 

Recently another attempt to describe the space-time evolution within the formation time was published in ref. [17]. 
The authors develop a quantum mechanical model to describe the expansion of the small wave packages. Their initial 
condition is motivated by the NRQGD approach from ref. [18]. The charmonium states are described as superposition 
of six charmonium states (4 5- waves and 2 P- waves). Diagonal and nondiagonal transitions due to collisions with 
nucleons arc taken into account. In contrast to the GVDM in the model of ref. [17] partons interact with the nuclear 
target within the formation time. Because this model also leads to a stronger suppression at smaller Feynman x, the 
cc pairs within the formation time have a smaller absorption cross section than the fully formed charmonium states. 
This is, because the formation time decrease with the decrease of the Lorentz factor of the cc pairs relative to the 
nuclear target. 

However, we restrict ourselves to the simpler model of refs. [14,2] because of simplicity and better control over the 
impact. For a quantum mechanical description like in ref. [17] a complete set of states is needed. The contribution of 
higher S and P states might be small, but the contribution of the continuum, given by open charm states, is unclear 
(remember that the t// is close to the DD threshold, and a mixing with DD may be relevant for the properties of 
large size configurations in the tp', cf. [19]). 

The NRQGD approach for the production of charmonium states is only valid if a large scale exists in the process, 
e.g. a high transverse momentum. For pt integrated data, where small transverse momenta dominate, this description 
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doesn't work. This is because [18] colored states arc produced that become color neutral by radiation of gluons. 
Without a large scale in the process only soft gluons can be emitted, which cannot transport quantum numbers (see 

e. g. [20]). Another problem in rcf. [17] is that a quantum mechanical model is unable to describe the emission of 
gluons. Therefore in this model the effect of the gluon emission on the cc-distribution is put in the initial distribution 
at the production point, though the emission of soft gluons takes a relatively long time. 

In our ansatz these deviations from quantum mechanics predicted by quantum field theories are taken into account 
in an effective way. In the expansion of cc pairs suggested in ref. [14,2] the cross sections of wave packages produced 
as small size objects increase with time, because the area of the color distribution increases. It is here not important 
if this increase is due to the motion of the partons or due to radiation of gluons. This scenario can be described like 
the diffusion of a QQ pair in statistical mechanics, because the total cross section of such a pair doesn't change due 
to gluon radiation if one sums over all channels, i.e., over the number of emitted gluons [21]. 

Ref. [22] claimed that the FNAL data on the J /tp photoproduction contradict to the expansion of cc pair discussed 
in ref. [14]. But ref. [14] discussed the space-time evolution of cc in the inclusive photoproduction of charmonium 
states (7 + A — > + X) for the energies where the coherence length is Ic ^ 2Ej/M^ <^ Ra- In the energy range 
of FNAL data the coherence length is comparable or exceeds nuclear radius Ra- At high energies the J/ip is not 
produced locally as it is the case at the SLAC energies considered in [14] and hence the analysis of [14] is not literally 
applicable - one needs to take into account both the formation of the cc before the nucleus and the increase of the cross 
section of cc-N interaction with energy. Besides ref. [22] discussed the quasielastic photoproduction of charmonium 
states {'J + A + p + {A — 1)). The selection of one specific channel means that there is no summation over 
all radiated gluons as for the inclusive process. Therefore the model of rcf. [14] should be modified to account for 
the suppression of gluon radiation, the expansion of the cc pair in the kinematics of FNAL photoproduction data is 
described better by the generalized gluon distribution, i.e., by the well understood QCD evolution ref. [23,24]. In this 
paper we discuss the total charmonium-nucleon cross sections at moderately large energies. Therefore the expansion 
model of ref. [14] should be applicable. 

C. Vector Dominance Model 

Somewhat different values for the cross section of the J/ip and V-''-nucleon interaction arise in charmonium models 
when the cross section is assumed to be determined by the radius of the color distribution within a hadron (tab. I) 
and also in the analysis of photoproduction data within the Generalized Vector Dominance Model (GVDM) [25] . The 
X-nucleon interaction can not be investigated with this model, because it has different quantum numbers than the 
photon. Therefore, parametrizations like in sect. II D are needed to make an educated guess for the x-nucleon cross 
section. 

The VDM takes into account only the direct diffractive production of the J/tp and the tp', while the GVDM accounts 
also for the non-diagonal transitions {tp' + N ^ J/ip + N and J/tp + N ^ tp' + N). The later are needed, because in 

photoproduction the particles arc produced as point like configurations and develop then to their average size. In a 
hadronic model like the GVDM this is taken into account in form of the interference due to the nondiagonal matrix 
elements. 

In the GVDM the photoproduction amplitudes /^^ and f^^' for the J/tp and the tp' are given by [13] 

e e 

f. ^ and f^i are the J /tjj — 7 and the tp' — j coupling and fvv are the amplitudes for the processes V + N ^ V' + N, 
where V and V are the J/tp and the tp' respectively. In the VDM the non-diagonal amplitudes with V ^ V are 
neglected. The importance of the nondiagonal transitions is evident, because the left hand side of eq. (7) is small. 
If it is neglected as a first approximation [13], then f^/^ = — j^/^,/, « 1.7 • /^^. And due to the CPT-theorem 
/jZ-'j/i = fjpjp'- The right hand side of eq. (7) is illustrated in fig. 1. 
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FIG. 1. Schematic illustration of the right hand side of eq. (7). 



To get a real description of this expansion, it would be necessary to have a complete set of hadron states, i.e., 
infinitely many bound S-wavc charmonium states plus the continuous spectrum (only S-wavcs have the same quantum 
numbers as the photon). For illustrating the physics, we restrict ourselves to the contribution of two states only. 

In ref. [25] the GVDM yields 8 10 mb for the ■^'-nucleon interaction cross section at SPS-energies. The sign of the 
photoproduction amplitudes f-y-^ of the J/ip and f^^^' of the ip'is positive relative to f-^i within the conventions 
of nonrclativistic charmonium models. The sign of the wave functions in the origin arc usually chosen to be positive. 
It is easy to see that a change of this convention would not change eq. (7). The J/i/i-nucleon interaction cross 
section at SPS-energies approximately 3.5 - 4 mb is used as input into the analysis of ref. [25]. The accuracy of such 
GVDM in predicting the tp'^ cross sections is not clear. The above calculation demonstrates that implementing color 
transparency leads to significantly larger cross sections of the tjj' N interaction. 

D. Charmonium-Nucleon Cross Section within the Chctrmonium models 

One possibility to evaluate the predominantly nonperturbative QCD contribution is to use an interpolation formula 
for the dependence of the cross section on the transverse size 6 of a quark-gluon configuration within the constituent 
quark model, 

(Tabs = c-b'^ (8) 

where b is the distance between the two constituent quarks, transverse to the collision direction. Such a form for the 
interpolation formula is also supported by pQCD, where the relation 

a{b) = ^b^a,{Q^)-xG{x,Q^) , (9) 

— 9/b^, was found [4] for small 6^ only. Beyond small b^ eq. (8) has no justification - it is merely an educated 
guess. 

The constant c can be adjusted within the constituent quark model: 6j = |rj and is the square of the pion 
radius. This radius is known from the Vector Meson Dominance Model to be r^r = ^ ~ 0.65 fm [26], nip « 770 MeV 
is the mass of the p-meson. This result is confirmed by measurements of the electromagnetic form factor [27] of the 
pion. Thus within this model 

c = a„N/bl = 25mb/1.06fm^ = 23.5mb/fm2 (10) 

We use in the following tr^jv = 25 mb. This is in good agreement with the data [28] for both energies discussed here 

Eiab{CERN) = 450 GcV and Eiab{Hera B) = 920 GeV, because the energy dependence of this cross section is not 
so strong in this kinematical region, b^ is increasing with energy as the cross section. Within soliton type models for 
the pion, the relation between b and the radius of the pion will be different. Therefore eq. (8) is model depended. 
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The X — N cross sections is calculated via: 



a = J • |*(r)|2dV . (11) 

where ^'(r) is the charmonium wave function. In our calculations we use the wave functions from two non-relativistic 
charmonium models: 1. A Cornell confining potential, 

52 r 

= + (2.34GeV-^)^ ' (12) 

sec [6] and refs. therein, and 2. a logarithmic potential 

Viogir) = -0.6635 GcV + (0.733 GcV) log(r • 1 GeV) (13) 

from ref. [7]. We chose two different charmonium models to exhibit the theoretical uncertainties of this approach. 
The resulting cross section are given in table I. 

This also resolves the puzzle why in photon- nucleus collisions a smaller (3-^4 mb) cross section for the J/ tp was 
fomid than in proton-nucleus collisions (6 7.5 mb): In pA collision 40% of the genuine J/ip^s are from decays of x 
mesons, which are strongly suppressed in 7A collisions, because the x has angular momentum L = 1, while J/ip and 
tp' has the same quantum numbers as the photon. 

TABLE I. The average square of the transverse distances of the charmonium states and the total quarkoniuni-nucleon 
cross sections a for two different charmonium models. For the x two values arise, due to the spin dependent wave functions 
{Im = 10, 11). The (Thard values are calculated with the eqs. (9) and (11). b < 0.35 is an upper limit for the integral in eq. (11). 
For the row a{hard,aUb) we used a parametrization like eq. (9) above 6 = 0.9 fm, because in this region Bjorken-a; becomes 
larger than 1 and the gluon distribution of eq. (9) cannot be defined in this region. 



cc/bb-state 


3/^ 


*' 


XclO 


Xcll 


<b' > (fm") 


0.094 


0.385 


0.147 


0.293 


a{Cornell) (mb) 


2.2 


9.1 


3.5 


6.9 


cr{log.) (mb) 


2.7 


13 


4.4 


8.7 


a(hard,b< 0.35 fm) 


0.42 


0.13 


0.41 


0.28 


a{hard, all b) 


2.2 


8.2 


3.7 


7.4 



E. Hard cross section 



We also calculate (t{X — N) assuming that pQCD is applicable. In this calculation we ignore the differences between 
bare quarks of the QCD Lagrangean and the constituent quarks, because these are nonperturbative QCD effects. The 
numerical calculations shown in Tab. I yields at CERN energies (7{X — N) values which are significantly smaller than 
in [5]. The ahard values [29] in Tab. I are also calculated with eq. (11), but here we integrate only over the region 
fm < 6 < 0.35 fm, because increases with b. b is the transverse distance between the heavy quarks transverse 
to the momentum of the heavy quarkonium. For the row a{hard, all b) we used a parametrization like eq. (9) above 
b = 0.9 fm, because in this region Bjorken-x becomes larger than I: the gluon distribution of eq. (9) cannot be defined 
in this region! For the calculation of the hard cross section, the a{b) in eq. (11) is given by eq. (9). We employ the 
CTEQ5L parameterizations of the gluon distribution functions in the proton [30] . 

The Bjorken x needed for the calculation of <Thard{X — N) is calculated by a; = / {2mMv), where — 9/6^, 
rriAr is the nucleon mass and v is the energy of the heavy Quarkonium state X in the rest frame of the nucleon. The 
calculation in Tab. I is done for a state X produced at midrapidity for SPS energies, but in the target fragmentation 
region for RHIC and LHC. One can see that the hard contributions to the cross sections are just a correction at SPS 
energies, but at RHIC energies both contributions become compatible, while at LHC the hard contributions dominate 
(we neglect here that the DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) equation might be violated [31]). 

The extrapolation of eq. (9) yields similar values as a parametrization like eq. (8). This is because the factor 
as ■ xG{x,Q^) depends only weakly on x and Q^, even in the region where <C f GeV^. This is no prove that 
parametrizations like eq. (8) can be applied to charmonium states, but at least it shows that there is no discontinuity 
between the pQCD result and the parametrizations for the soft regime. 
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III. EXPERIMENTAL RESULTS 



A. Comparisson with CERN Data 

In the lefthand side of fig. 2 we show a comparison between calculations with different cross sections and different 
expansion times and the NA50 data [1] for pA collisions and the NA51 data [32] for pp and pD collisions for the cross 
section of -0' N interaction vs. the mass of the target. The y-axis shows iJ^^cr^' /A where i?^^ is the branching ratio 
for the decay of the ip' into dimuon pairs, and a^i is the production cross section. The "5.1 mb, instant formation" 
curve in fig. 2 (Ihs) is the fit of the NA50 collaboration to their data. Instant formation means that they assumed 
that the tp' is produced with the full cross section and not as a point like particle as in the description of this paper. 
(Note the the NA50 collaboration fitted B^j^f^a^' /any, where DY means Drell-Yan, therefore we multiplied there fit 
with the Drell-Yan cross section in pp collisions measured by NA51). 

The "8 mb, tf = 0.45 fm" curve is the eye-ball fit of the model described in this paper. For the comparison with 
the data we need the production cross section of the ip' in pp collisions as input. We used here the average of the pp 
and pD data of the NA51 collaboration. The value of 8 mb agrees well with the model parametrizations discussed 
in the sections II A and II B. However, we compare also with the calculation with the parameters of ref. [2], i.e., 
a{'tp'N) = 20 mb and tf = 0.6 fm. For this comparison we used the production cross section of the i/j' in pD collisions 
divided by two as input. This is also close to the value of the NA50 fit. One can see in the lefthandsidc of fig. 2 that 
the calculation with these parameters is also in good agreement with the data. A value for a{ip' N) of the size of 20 
mb is favored by the nucleus-nucleus data as shown in ref. [8] . 

In the righthandside of fig. 2 is plotted B ^^a{J / ip) / <y{DY) the NA50 data [33] for PbPb collisions vs. the transverse 
energy Et, a measure for the centrality of the collision. The result of the calculation within the UrQMD (Ultrarel- 
ativistic Quantum Molecular Dynamics) model [34] is also shown in form of the histogram in fig. 2 (rhs). In the 
calculation (for more details see ref. [8]) for the charmonium nucleon absorption cross sections the values from ref. [2] 
were used. In this calculation the charmonium states can interact also with secondary produced particles. Based on 
the additive quark model the cross section for charmonium baryon interactions was assumed to be the same as the 
charmonium nucleon cross sections and for charmonium meson interactions two thirds of the charmonium nucleon 
cross sections. 

The calculation for the J /ip agrees well with the data. The calculation for the ■0' underestimates the data. However 
it is not understood, if this is due to the high value of a{ip'N) = 20 mb, or if nondiagonal transitions like in sect. II C 
should be taken into account in AA collisions, too. 

The value of 8 mb for the ip' nucleon cross section from the fit to the pA is smaller than the theoretical estimate 
20 mb of ref. [2]. This is because in ref. [2] a formation time of 0.6 fm was chosen for the -0': while we used here 0.45 
fm, the radius of the ip' given by the charmonium models. The fact that the formation time is not known very well is 
another uncertainty. Further uncertainty comes from using the diffusion model of expansion at distances comparable 
to the scale of the soft interaction. Within the error bars the V'-nucleon cross section extracted from these pA data 
and the prediction of the GVDM, discussed in sect. II C, are qualitatively similar. However, further data are needed 
to learn more about this cross section. 



B. Predictions for HERA B 

Fig. 3 shows our predictions of S{J/'ip) lefthand side (Ihs) and S{ip') (rhs) for HERA B vs. Feynman x for pA 
collisions at Eiab = 920 GeV. Two different nuclear targets. Carbon {A = 12) and Tungsten with [A = 184), are 
shown for the two different sets of cross sections resulting from the different charmonium models. 

We used here the cross sections of the constituent quark model shown in tab. I in sect. II D, because this model is the 
only one which predicts also nucleon cross sections. We calculated the suppression factor S for the ip' additionally 
with a cross section of 20 mb and an formation time of 0.6 fm to estimate theoretical uncertainties. These are the 
values discussed in ref. [2]. 
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FIG. 2. On the (Ihs) are plotted the Bfj,^a^i /A values extracted from calculations with different absorption cross sections 
and different formation times, the NA51 data for pp and pD, and the NA50 data for pBe, pCu, pAg, and pW vs. the mass A 
of the target. And the (rhs) shows Bi_i,fj,a{J /ip) / a{DY) with the absorption cross sections from ref. [2] and the NA50 data for 
PbPb collisions vs. the transverse energy Et- 
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FIG. 3. SiJ/tp) (Ihs) and SO^') (rhs) are shown vs. xp for Carbon and Tungsten at HERA B energies {Eiab = 920 GeV). 

Note that we plotted in fig. 3 the result of our model up to a Feynman x of 0.2 to show the behavior of the model. 
However, it is known from E866 [3] and NA3 [35] data that there are additional suppression mechanisms at positive 
Feynman x, e.g. parton energy loss and nuclear shadowing. 
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IV. CONCLUSIONS 



The new data of the NA50-collaboration [1] and the data of the E866-collaboration [3] prove that the ■!/''-iiucleon 
cross section is much larger than the J/ip-miclcon cross section. This is in agreement with the photoproduction data 
for these charmonium states as discussed in the framework of the GVDM in section II C. This confirms the QCD 
prediction that the strength of hadron-hadron interactions depends on the volume occupied by color. 

Within the assumption that charmonium states are produced as point like white states, we demonstrated that the 
data [1] can be fitted with a ip'-micleon cross section of (T{tp'N) « 8 mb. We discuss different models for this cross 
section and show that they prefer a -0'-nucleon cross section of (7{tp'N, model) = 9 13 mb. However, due to the 
large experimental errors we conclude that the data and the QCD-motivated models agree, further data with higher 
accuracy and covering a larger rapidity range are needed. 

This cross section will be measured soon in proton-nucleus collisions at HERA B at an energy of Eiab = 920 GeV. 
The advantage of this experiment is that it covers a larger range of Feynman-XF, especially in the negative xp region. 
In this region effects due to the formation time of the hadron will be less important. These data will give new 
information about the relations between the size of mesons and the strength of hadron-hadron interactions as well as 
the possibility to test the existing models for this relation. 
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